Astronomy <fe Astrophysics manuscript no. H4034 


February 2, 2008 


(DOI: will be inserted by hand later) 





Ionization structure in the winds of B[e] supergiants 
I. Ionization equilibrium calculations in a H plus He wind 
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Abstract. 

The non-spherically symmetric winds of B[e] supergiants are investigated. An empirical density distribution is 
chosen that accounts for the density concentrations and ratios derived from observations, and our model winds 
are assumed to contain only hydrogen and helium. We first calculate the approximate ionization radii for H and 
He and compare the results with the ionization fractions calculated from the more accurate ionization balance 
equations. We find that winds with a r~ 2 density distribution turn out to reach a constant ionization fraction as 
long as the wind density is low, i.e. in polar direction. For the high density equatorial regions, however, we find 
that the winds become neutral just above the stellar surface of the hot and massive B[e] supergiants forming a 
disk-like neutral region. In such a disk molecules and dust can form even very near the hot central star. 

Key words, supergiants - Stars: mass-loss - Stars: winds, outflows - Hn regions 



1. Introduction 

B[e] supergiants are massive and luminous post-main se- 
quence stars with strong non-spherically symmetric winds 
induced probably due to the rapid rotation of the stars. 
Observations of dust and molecules indicate that there 
must exist some neutral regions near these stars that are 
shielded from the ionizing and dissociating stellar radi- 
ation field. The goal of this paper is to show that in the 
equatorial regions where the slow, high density wind forms 
a disk-like structure, hydrogen becomes neutral just above 
the stellar surface and can therefore perfectly shield the 
disk material to allow molecules and dust to be formed 
near the hot central star. 

Photoionization of gas around hot stars is since long 
a well known process leading to the beautiful observable 
emission objects like Hn regions and planetary nebulae. 
The central stars thereby ionize the material that has 
been blown away by the star itself in a more or less dense 
wind. The sizes of such ionization regions, also called the 
Stromgren spheres, can be calculated by balancing the 
number of ionizations by the number of recombinations 
and usually a very sharp outer edge of the ionized region 
is found. These calculations are normally performed for 
nebulae with a constant electron temperature and den- 
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sity (for an overview see Osterbrock ( 1989) and references 
therein) . 

O and B stars have a stellar wind which for a spher- 
ically symmetric and steady radial outflow leads to a 
v(r) r wind density distribution. So the calculation 
of the ionization structure in these winds is more com- 
plicated than for a homogeneous cloud. In an early work 
by Drew (1985) a radiative ionization and thermal equilib- 
rium model for the wind from P Cygni was calculated, and 
later on e.g. Pauldrach (|1987f> . Drew ([19891, MacFarlane 
et al. l|1994f) . and Pauldrach et al. (|1994fl calculated in 
great detail the ionization structure in the spherically sym- 
metric winds of hot luminous stars using the ionization 
balance equations. In addition, they solved the thermal 
equilibrium equations to calculate the electron tempera- 
ture in the wind. For this, line cooling by collisionaly ex- 
cited heavy elements is the most important cooling mecha- 
nism and the authors took a huge amount of different lines 
into account which results in an enormous set of equations 
that has to be solved. 

In this paper we calculate the ionization structure in 
the non-spherically symmetric winds of B[e] supergiants. 
For simplicity, these winds are supposed to contain only 
hydrogen and helium. We especially focus our investiga- 
tions on the equatorial regions where we expect to find 
neutral material near the central star. In our calculations 
we make use of several approximations which result in an 
upper limit for the distance of the equatorial neutral re- 
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gion, because we want to show that efficient shielding of 
dust at larger distances can occur. 

2. The winds of B[e] supergiants 

B[e] supergiants are thought to be post main-sequence 
stars but their real location within a clear evolutionary se- 
quence is still under debate (see e.g. Schulte-Ladbeck ri997l 
Lamers et al. 120031 . Up to now about 15 B[e] supergiants 
are known in the Magellanic Clouds, and a few candidates 
are known in the Milky Way although the classification as 
a B[e] supergiant in the Milky Way is very difficult due 
to the lack of good distance and therefore luminosity de- 
terminations (Lamers et al. I1998fl . The situation is more 
favorable for B[e] supergiants in the Magellanic Clouds for 
which the distances are known. 

The B[e] supergiants in the Magellanic Clouds have 
been observed in great detail and many spectra in dif- 
ferent wavelength regions are nowadays available. Most 
interesting, these stars were found to have a so-called hy- 
brid spectrum (Zickgraf et al.ll985t which means that on 
the one hand the UV resonance lines of highly ionized 
atoms show broad P Cygni profiles indicating a wind of 
low density and high velocity, while, on the other hand, 
forbidden and permitted emission lines of low ionized ions 
and of neutral atoms arising in the optical spectrum indi- 
cate a wind of high density and low velocity. The winds 
of B[e] supergiants must therefore consist of two different 
components and Zickgraf et al. l|1985fl proposed a picture 
of these winds containing a hot and fast radiation-driven 
CAK-type wind (Castor et al. H975[) in polar direction and 
a much cooler, slow and dense wind in equatorial direc- 
tion which might form a disk-like structure. The density 
contrast between equator and pole is of order 10 2 — 10 3 , 
and from the statistical viewpoint a mean opening an- 
gle of 20° — 40° of the disk has been derived (Zickgraf 
1990] The reason for the non-spherically symmet- 

ric winds is probably the rapid rotation of these stars. 
Rotational velocities in the equatorial region of the star 
of about 40 — 70 % of the break-up velocity have been de- 
rived from their spectra (Zickgraf et al. I199fi|> . Pelupessy 
et al. H2000|) and Bjorkman (1998) have argued that the 
disk is probably due to a combination of rotation induced 
bi-stability and rotation induced wind compression (see 
Lamers & Casinelli ITM?! Chapter 11). 

Other observations confirm the existence of a neu- 
tral disk. These are emission from molecules like CO 
(McGregor et al. I1988al I1988b|) and TiO (Zickgraf ct al. 
1989) and emission from cool and warm dust with strong 
contributions to the spectral energy distribution (Zickgraf 
et a,l. [ll^[TT^m^rTM2l Gummersbach et al.[TH23). 
Also some typical dust emission features have been ob- 
served (e.g. Waters et al. Allen HWR Zickgraf et 
al. I1986|l . The dust component in these disks is also con- 
firmed via polarization measurements (Magalhaes 



Detailed observations of the B[e] supergiants with different 
inclination angles show that the polarization is strongest 
for the edge-on systems. In addition, there is a correla- 




100 



150 



'i 1 1 1 1 1 




'I 1 1 1 I 1 1 ■ I 1 1 1 I 1 1 1 



2 4 6 8 100 0.2 0.4 0.6 0.8 1 1 
F m (e)/F m . pol v.(S)/v., pol 



400 800 
Nh(«)/Nh.poi 



Fig. 1. Angular distribution of the mass flux, velocity and 
hydrogen density on the stellar surface. The parameters 
are b\ — 1; 62 = —2; s — 100 (dashed lines); b\ — 0.5; 62 = 
-1.5; s = 100 (dotted lines); and 61 = 1; b 2 = -2; s = 10 
(solid lines). 

tion between the polarization and the [H — L] color of the 
stars, indicating that the observed infrared excess emission 
is due to circumstellar dust concentrated in the equatorial 
region. This dust disk is the origin of the major part of 
the observed intrinsic polarization. 

3. The wind model 

For a star with a stationary spherically symmetric wind, 
the mass loss rate is related to the density and the ve- 
locity at any point in the wind via the equation of mass 
continuity 



M = 47rr 2 p(r) v(r) 



(1) 



The mass flux per unit surface area from a spherical star 
is 



F m = 



M 
4vri? 2 



p(R*)v(R*) 



(2) 



For a non-spherically symmetric wind the mass flux is a 
function of latitude, F m — F m (9), and the mass loss rate 
is 



it/2 

M = 4ttRI [ F m (9)sm(6)d6 



(3) 



where we assumed that the star is a sphere. The hydrogen 
density distribution in such a wind is 



1 fR*YF m (0) 



limn V r ) v(9,r) 
which we parametrize in the following form 



1 f R*\ F, 



pmu \ r J v po i(r) 



m,pol 



^gb-(sin0) s 



(4) 



(5) 



The term 10 & with 6 = 2 — 3 reflects the density contrast 
between the equator and the pole, and the term (sm6*) s 
describes the stronger density concentration with increas- 
ing s towards the equator. The parameters F m!po i and 
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Table 1. Model wind parameters. Listed are the electron temperature, T e , mass fluxes in polar and equatorial direction, 
-Fm.poi and -F mjCq , mass loss rate of the star, M, terminal velocity in equatorial direction, foo^q, and the resulting 
equatorial hydrogen density on the stellar surface, -/Vjj.eqC-R*)- The velocity in polar direction is assumed to be in all 
cases 2000 km/s which is a typical value in a CAK-type wind. 



Model 




F u 


,pol 


-^m,eq 


M 


^oo ,eq 


b 


s 


iV H , e q(-R*) 




[K] 


[g s" 1 


cm" 2 ] 


[g s _1 cm -2 ] 


[M /yr] 


[km/s] 






[cm" 3 ] 


A 


1.0 ■ 10 4 


1.5 • 


io- 7 


1.5 • IO -0 


2.4 • 10 _b 


20 


3.0 


100 


3.3- 10 11 


B 


1.5 • 10 4 


1.5 ■ 


io- 7 


1.5 ■ 10~ 6 


2.4 • IO -6 


20 


3.0 


100 


3.3- 10 11 


c 


1.0 ■ 10 4 


1.5 • 


1Q -6 


1.5 ■ 10~ 5 


2.4 • 10~ 5 


20 


3.0 


100 


3.3- 10 12 


D 


1.0 ■ 10 4 


1.5 • 


IO -8 


1.5 ■ 10" 7 


2.4 ■ 10~ 7 


20 


3.0 


100 


3.3 ■ 10 10 


E 


1.0 ■ 10 4 


1.5 • 


io- 7 


4.9 ■ 10~ 7 


1.8 • 10~ 6 


63 


2.0 


100 


3.3 ■ 10 10 


F 


1.0 • 10 4 


1.5 • 


IO" 7 


1.5 ■ 10" 6 


4.0 ■ 10~ 6 


20 


3.0 


10 


3.3- 10 11 



Vpoi(r) are the mass flux and wind velocity in polar di- 
rection which can be derived from the observations. Their 
latitude dependence is taken to 



^ m (0)=^m,pol-lO MSine) 

and 

v(6,r) =v pol (r) -10^^ 



(6) 



(7) 



with b — b\ — 62 • This implies that we have adopted the 
same shape of the velocity law, v(r)/voo, for all latitudes. 
In Fig.n]we show the variation of the mass flux, velocity 
and hydrogen density along the stellar surface for different 
values of 61, 62, and s. 

The radial velocity dependence in a stellar wind is nor- 
mally described by the /3-law 



v(r) 



r 







(8) 



where Voo is the terminal velocity. For hot stars, /3 is typ- 
ically of order 0.8, and the wind reaches its terminal ve- 
locity at a distance of a few stellar radii. 

For our calculations we make the following assump- 
tions, whose validity and influence on the results are dis- 
cussed in Sect. EJ 

(i) we set the electron temperature in the wind constant, 

(ii) we take the velocity in the wind to be independent of 
distance, i.e. v(6,r) = v oo (0). 

For this paper we choose a star with a mass of ap- 
proximately Af» = 30 Mq, a radius of i?* = 82 i?Q, and 
a radiation temperature of the photosphere in the Lyman 
continuum of T rad = 1.7 • 10 4 K which coincides with a 
Kurucz model atmosphere fKurucz ll979f) for a star with 
an effective temperature T c g ~ 22 500 K and logg = 3. 
The luminosity of this star is = 1.5 • 10 6 L Q which cor- 
responds to the most luminous B[e] supergiants observed. 

For the wind we calculate different models whose pa- 
rameters are summarized in Tablc^ 

4. Ionization radii in a wind with non-spherical 
density distribution 

In this section we discuss the calculation of ionization 
radii. We adopt a similar procedure as for the classical 



Stromgren sphere calculation but use the non-spherical 
density distribution given by Eq. (JSJ) ■ We first derive the 
general formulae. In Sect. 14. 21 we discuss the ionization 
radii in the optically thin case while in Sect. l4.3l the on- 
the-spot approximation is used. Both models arc with the 
assumption that the ionization is complete within the ion- 
ization radii (as usually adopted). This simplification al- 
lows us to get a first impression of the overall shape of the 
ionized region around a B[e] supergiant fSect. l4~^|l . The 
physically more relevant case of partial ionization is dis- 
cussed in Sect.|21 

4.1. General formulae 

We restrict the calculations to a wind containing hydrogen 
and helium. Since the following equations have the same 
shape for all angles 8 we skip the angular dependence of 
the densities and write down the equations only as func- 
tions of the radial distance. In addition, we assume that 
the fraction of He ++ is negligible. 

The ionization radius is calculated by balancing the 
number of ionizations which depends an both, the stellar, 
Jl, and the diffuse radiation field, , and the number of 
recombinations per cm 3 and per s. For H and He we can 
therefore write 



iV H o(r) J a„(H° 



hv hv 
N m (r)N c (r)a A (B°,T) 



dv 



00 

AWM / ^ (He ) 



4ttJ, s , AirJ?, 



hv hv 
N Hc+ (r)N c (r)a A (Re°,T) 



dv 



(9) 



(10) 



with N as the number density of the individual ions and 
neutral atoms, whereby 7Vh c = y ■ and y is the He/H 
ratio by number which we set to 0.1. The terms O!a(H , T) 
and aA(He°,T) are the radiative recombination coeffi- 
cients of hydrogen and helium for recombinations into all 
levels (see Fig.|2J, a„(H°) and a„(He°) are the ionization 
coefficients, vq and v^ are the ionization frequencies of hy- 
drogen and helium, and the electron density, N c , is given 
by N e (r) = N B+ (r)+N Bc+ (r). 
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Fig. 2. Hydrogen and helium recombination coefficients. 
The dots for H refer to values given by Seaton i|1959|l , the 
dots for He are taken from Osterbrock Ijl989|l . and the 
lines give fits according to interpolation calculations. 



Taking into account the absorption of stellar contin- 
uum photons along the direction to the star, the mean 
intensity of the stellar radiation is 



4irJl = 4nW(r)B v (T iad )e- T " 



(11) 



where W(r) = |(1 — y/l — (R*/r) 2 ) is the geometrical 
dilution factor, T ra d is the radiation temperature of the 
photosphere, and t v is the radial optical depth defined as 



T v (r)= I a u (E°)N n o(r')dr' 
for i/Q < v < v-2, and 



(12) 



r„(r) = / (a,(H°)^ H «(r') +a,(He°)iVHeo(r')) dr' (13) 

for V2 < v. We know that Eq. is only correct as long 
as the star can be considered as a point source. But for 
our purpose this approximation is also good enough for 
distances near the star (see Sect.[5J. 

The diffuse radiation field is much more complicated to 
calculate because at every point in the wind one must in- 
tegrate the diffuse radiation over the complete wind zone. 
To do so, one needs the information of the structure of 
the wind zone, i.e. information about where in the non- 
sphcrically symmetric wind exist neutral regions, what is 
the absorptivity and emissivity of these neutral regions, 
etc. Therefore, we restrict our investigations to two ex- 
tremes: (a) the optically thin case, where all ionizing ra- 
diation produced by recombinations leaves the wind zone 
without being absorbed by a neutral atom, and (b) the op- 
tically thick case, better known as the on-the-spot approx- 
imation, which states that every ionizing photon produced 
in the wind via recombination is immediately absorbed by 
a neutral atom nearby, i.e. on-the spot. In reality the mean 
intensity will be in between the values predicted for these 
two extremes. 



4.2. Non-spherical ionization radii in the optically thin 
approximation 

In this case, the diffuse radiation field is zero, and the 
ionization equilibrium equations for H and He, Eqs. 
and reduce to 



47iW(r)(l - <te(r)) 



ff,(T rad K(H°) 
hv 



M dv 



= fe(r) N B (r) (q B (r) + yqSc (r)) «a(H°, T) (14) 



and 



47rW(r)(l- te .(r)) / B *( T »dK(He ) e -Mr) dv 
J hv 

- <ZHcM N H (r) (qn(r) + yqn c {r)) a A (He°, T) (15) 

where we replaced the densities by density fractions 
N n +(r) = qn(r) • Nn(r) &ndN n o(r) = (l-«H(r))-JV H (r), 
and the same for He. 

To calculate the distance from the star at which the 
gas turns from ionized into neutral, i.e. the ionization ra- 
dius R s (0), we assume (as usually adopted in classical 
Stromgren sphere calculations) that all atoms are ionized 
whithin the ionization radius, and all atoms are neutral 
outside, which means that 

/ 1 5 r < Rf , / 1 ; r < Rf c , . 

9h = i o . r > r h and fflc = <j n _ ^ D Hc (16) 



; r > Ri 



Integrating now both sides of Eqs. I)14|l and (|15J) over the 
radius, and taking into account the different sizes of the 
ionized regions for H and He with i?f c < Rf, the ioniza- 
tion balance equations become 



oo oc 



4tt JW{r) 



hv 



W dvdr 



a A (H°,T) 



R" 



R, 



Ns(r)(l +y)dr+ J Nn(r) dr 



(17) 



for hydrogen, and 



oo oo 



_RHc V2 



a A (He°,T) J N H (r){l+y) dr 
R, 



(18) 



for helium. The right-hand side of Eq. (|17|) splits into two 
integrals according to the two different regions that have 
to be considered: (i) the region between the stellar surface 
and the He ionization radius where both, H and He are 
singly ionized, and (ii) the region between the He and H 
ionization radii where only H is ionized. Obviously, the 
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corresponding integration limits for the optical depth cal- 
culation must be adapted correctly. From Eqs. ljl"7|l and 
<|18|) the ionization radii for H, R^, and He, R^°, can be 
calculated for every angle 6. 

These equations look quite different from the classical 
Stromgren sphere formula. But since they are more gen- 
eral, the classical Stromgren sphere, as a special solution, 
can be found easily as shown in Appendix A. 

4.3. Non-spherical ionization radii in the on-the-spot 
approximation 

The on-the-spot approximation for optically thick winds 
assumes that every ionizing photon generated by recom- 
bination will be absorbed very near its generation loca- 
tion. The diffuse field can be split into several parts and 
we consider diffuse photons generated by (a) recombining 
hydrogen atoms into the ground level, (b) recombining he- 
lium atoms into excited levels, and (c) recombining helium 
atoms into the ground level. In this paper we take into ac- 
count neither ionization from excited levels nor collisional 
ionization and recombination. Both effects will strongly 
influence the ionization balance but we expect that they 
become really important only for stars with much higher 
radiation temperature in the Lyman continuum than the 
17 000 K we are dealing with. A full solution of the level 
population equations together with the ionization balance 
equations, where these effects are included, is presently 
under investigation and will be published in a forthcom- 
ing paper. 

While all of the diffuse photons generated by the three 
mentioned processes are able to ionize hydrogen, only pho- 
tons resulting from the third process can ionize helium. To 
complete the ionization equilibrium equations for H and 
He in the on-the-spot approximation, we have to add the 
term 

flk(r) (^hM+^HcM) [g H (r)ai(H°,r) 
+yqMr) (£ai(He ,T) +pa B (He°,T))] (19) 

that describes the additional ionizations of H due to the 
diffuse radiation field, to the left-hand side of Eq. fH| . and 
the term 

N H (r) (m(r) + yqnc(r)) </<?He(r) (1 - £) ai(He°,T) (20) 

for the additional ionization of He due to the diffuse 
radiation field, to the left-hand side of Eq. I|15|) . With 
the electron density in terms of the q- values, N e (r) = 
Nn(r) (q-ii{r) + yqn c (r)), the ionization balance equations 
that must now be solved can be written in the following 
form: 

oo 

ATrW(r) (1 - g H (r)) / ^( r ^K(H°) e -r„(r) dv 

+ N c (r) [g H Wai(H ,T) + y(? HcW (^i(He°,T) 
+ pa B (He°,T))] = N e (r) q H (r) a A (H°, T) (21) 



for hydrogen, and 

oo 

AnW(r) (1 - q He (r)) J ^Z^&l £ -Mr) dv + Nc(r) 

•y<ZHc(r)(l-Oai(He°,T) = ^ c (r) gH cWaA(He ,r) (22) 

for Helium. Here, q:b is the recombination coefficient to 
all levels except of the ground level, i.e. qb = a a — &i, 
and £ is the fraction of recombining helium to the ground 
level which is absorbed by H atoms. The remaining frac- 
tion 1 — £ is absorbed by He atoms, p is the fraction of 
photons generated by recombination of He to all excited 
levels that can ionize H. This fraction strongly depends 
on the electron density, because, if the electron density, 
N e (r), at the location r is high enough, then collisions of 
electrons with the excited He atoms will depopulate the 
levels and no or only a few photons are produced. The 
critical density, N c , for collisional deexcitation is defined 
as the density at which collisional transitions are equally 
probable with radiative transitions. The values for the ra- 
tio p are either p ~ 0.96 if N c (r) <C N c (r) or p ~ 0.66 if 
iVe(r) > N c {r) (see Osterbrock 11555^ . 

Making again the assumptions that all atoms are ion- 
ized within and neutral outside the ionization radius, ig- 
noring the absorption by H in the He + zone by setting 
£ = 0, and integrating the complete ionization equilib- 
rium equations over the radius yields again R s (0) for H 
and He. 

4.4. Results 

The ionization radii, R s (9), calculated for the different 
model parameters (see TableQJ, are presented in Fig.|3]for 
the optically thin approximation and in Fig.0|for the on- 
the-spot approximation. In addition, Figs.0 and show 
the ionization radii of the models D and E but for a much 
larger distance. 

The most important result is that for all of the models 
there is a disk-shaped neutral region around the star, and 
for the high-density models (A, B, C, and F) these neutral 
disks can even exist at or just above the surface of the hot 
star. 

Due to the higher number of ionizing photons avail- 
able in the on-the-spot approximation we would for this 
method expect the ionization radii to be larger compared 
with the optically thin results, because e.g. in the case of 
He the size of the ionization radius depends only on the 
recombination coefficient. Comparing Figs.|3] and 0] this 
effect is not obvious for the equatorial direction due to 
the choosen scaling. But the ionization radii calculated in 
the on-the-spot approximation are indeed larger as can be 
seen by comparing the results for a more polar direction 
or by comparing the Figs .|5] and |S| The real location of the 
ionization radius should lie in between the two boundaries 
defined by the optically thin and the optically thick ap- 
proximation. 
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Fig. 3. Ionization radii in equatorial regions of H (dashed 
lines) and He (solid lines) for the different model calcu- 
lations in the optically thin approximation. The stellar 
surface is shown as the dotted line. 



We take model A as the reference model and discuss 
now the effects of the change in the different free parame- 
ters on the size of the ionization radii. The higher electron 
temperature in model B leads to a slightly larger ion- 
ization radius because the recombination coefficient and 
therefore the number of recombinations decreases with in- 
creasing electron temperature (see Fig-EJ) . The influence of 
the electron temperature is, however, not very significant. 

The series of the models C, A, and D displays the 
behaviour of the ionization radii for a decreasing mass flux, 
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Fig. 4. As in Fig.[3|but for the on-the-spot approximation. 



i.e. a decreasing hydrogen density in equatorial direction. 
Consequently, the ionization radii of H and He increase 
drastically because (i) much less recombinations occur, 
and (ii) for a lower density, i.e. a lower optical depth in 
direction to the star, the ionizing photons can penetrate 
much deeper into the wind. While the density in model 
A is still high enough for effective recombination and for 
an effective absorption of the ionizing photons near the 
stellar surface, especially the H ionization radius in model 
D is shifted to a much larger distance (see Figs. [S] and [6J. 

Not only the change in mass flux leads to a change in 
density. Also a lower value of 6, the density ratio between 
equator and pole, results in a lower density and leads to 
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Fig. 5. Optically thin approximation for the models D and 
E but for a much larger distance. Shown are the ionization 
radii for H (dashed line) and He (solid line). 
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Fig. 6. On-the-spot approximation for the models D and 
E but for a much larger distance. Shown are the ionization 
radii for H (dashed line) and He (solid line). 

a shift of the ionization radius to a much larger distance 
(model E, see also Figs. El and EJ). 

And finally, the smaller value of s in model F, which 
means that the concentration of the density to equatorial 
directions is less strong, leads to a much broader neutral 
disk-like structure, as expected. 

5. Ionization structure of H and He 

In the previous section we calculated the ionization struc- 
ture using the approximation that the wind is fully ionized 
at r(0) < Rs(8) and neutral at r(9) > R s {9). This is a se- 
vere simplification because the region where the wind is 
partly ionized may be substantial. Therefore, in this sec- 
tion we calculate the ionization structure of the wind more 
exactly. As concluded from the Figs .0 and 0] the neutral 
region in equatorial direction can already be reached at 
or near the stellar surface, as long as the wind density is 
high enough. For these high density models, the on-the- 
spot approximation is a very good tracer for the H and 
He ionization radii and we use this approximation further 



on to calculate the ionization structure in these winds in 
equatorial directions. 

We start with the ionization balance equations in the 
on-the-spot approximation derived above, i.e. Eq. (|21|l for 
H and Eq. (|22|) for He. The optical depth is given by the 
Eqs. O and $F$. From Eqs. J2U and (T22JI the ionization 
fractions for H and He can be calculated radially outwards 
in the wind whereby the optical depth at each point in the 
wind between the stellar surface and the distance r must 
be known. The inclusion of the optical depth along the 
direction to the star into the calculations leads to results 
for the ionization radii which are more accurate than in 
the approximate calculation of the previous section. 

In Fig.Hwe plot the fractions of ionized and neutral 
helium calculated for the different models of TablefJ] The 
plots in the three panels are thereby grouped according to 
the values of the equatorial hydrogen density, i.e. model C 
in the top panel has the highest density, while the models 
D and E in the bottom panel have the lowest density. 
Also shown, as a thick vertical line, are the locations of 
the approximate ionization radii extracted from Fig.0J 

At first, in the top and middle panel the distance of the 
approximate ionization radius is much larger than for the 
more accurate ionization fraction calculations. The reason 
for this is the optical depth which is not really zero as 
assumed in the approximate ionization radius calculations. 
Instead, there exists always a very small fraction of neutral 
atoms which results in a non-negligible attenuation of the 
stellar radiation field. For example, in the highest density 
model C where the He neutral region nearly coincides with 
the stellar radius, the approximate method results in an 
ionization radius that is about a factor of 8 larger. This 
strong descrepency can be understood, if we look at the 
fraction of neutral He which turns out to be a few times 
10 -3 at the stellar surface. Such a large neutral fraction 
results in a substantial optical depth already at the stellar 
surface. 

For the models A, B, and F, i.e. the models with the 
intermediate density, the He neutral region starts at larger 
radius and the neutral fraction inside this region is much 
smaller (only about 10 -5 -10 -6 ). The difference between 
the approximate and the real ionization radii is therefore 
much smaller, about 25%. 

In the lower panel of Fig.[7|the results for the models 
with the lowest density, D and E, are shown. Here, the 
situation is different: the exact calculations lead to a larger 
distance of the neutral region than the ionization radius 
calculations. These seems not to be physically and shows 
us that in these two models the assumption of the validity 
of the on-the-spot approximation no longer holds because 
the density in the models D and E is already too low to 
guarantee that every photon produced by recombinations 
will be absorbed in the wind. 

The main result of this comparison in Fig.0is that the 
approximate calculations, by neglecting the optical depth 
effects inside the ionized regions, overestimate the real lo- 
cations of the ionization radii. The same behaviour can 
also be seen in Fig-IHl which displays the fractions of ion- 
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Fig. 7. He ionization structure in equatorial direction 
(9 = 90°) for the different models. The models are grouped 
according to their density in equatorial direction: the up- 
per panel belongs to the highest density, the lower panel 
to the lowest density models. The shift of the transition 
radius to higher distances for lower densities is clearly vis- 
ible. The solid lines represent qn e + and the dashed lines 
Que — 1 — <ZHc+ ■ The thick vertical line is the location of 
the ionization radius from Fig.^J 



ized and neutral hydrogen together with the locations of 
the approximate ionization radii. 

Interestingly, when calculating the ionization fractions 
for the models D and E we could not find the transition 
from ionized to neutral H, but instead the calculations re- 
sulted in a constant ionization fraction contrary to what 
was found for the approximate calculations. The same re- 
sult was achieved for calculations in polar directions. This 
discrepancy can easily be understood by inspecting the 
balance equations l|21l) and (|22() : If we summarize on the 
right hand side the terms that contain the electron density, 
the left hand side is proportional to W(r)-Nn(r) while the 
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Fig. 8. As Figure [7| but for the H ionization structure. 
Models D and E are not shown because due to their low 
densities no transition from ionized to neutral occurs; the 
ionization fraction in these winds becomes constant. 

right hand side is proportional to N^(r) ~ r~ 4 . Since for 
large distances, the dilution factor becomes proportional 
to r~ 2 , the left hand side also becomes proportional to 
r -4 which then drops out of the equations. The only r de- 
pendent term that is still in is the optical depth. But for 
large distances with only slowly varying q values the opti- 
cal depth becomes constant. For large distances from the 
star, the ionization balance equations therefore become in- 
dependent of r, and the ionization fractions become con- 
stant. As long as there is no change in the density distribu- 
tion, e.g. when the wind reaches the shockfront where the 
wind material hits the interstellar medium, there will be 
no neutral region in the polar directions of the wind but 
all material stays ionized with a constant ionization frac- 
tion. In the approximate calculation of the ionization radii 
we make the assumption that the neutral region starts at 
a distance smaller than infinity, but this assumption does 
not hold for low density winds with a r~ 2 density distri- 
bution. Nevertheless, in winds with a constant density, as 
might be the case in planetary nebulae, the assumption is 
correct because in a constant density wind the ionization 
balance equations never become independent of r. 



6. Discussion 

The motivation for this paper was to find the location of 
neutral material in the equatorial regions around hot and 
massive stars. In our model calculations we make some 
assumptions that seem not to be valid for the regions of 
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our interest, i.e. very near the stellar surface. In this sec- 
tion we therefore discuss these assumptions and show that 
they are indeed justifiable also near the stellar surface. 

1. The constant electron temperature. In a test cal- 
culation we computed with the help of the thermal 
equilibrium equations for a pure H plus He wind the 
electron temperature distribution in equatorial direc- 
tion. We found that on this short distance between 
stellar surface and neutral region the electron temper- 
ature varies only very slightly, so that our assumption 
of a constant electron temperature seems to be valid, 
at least for the equatorial regions. 

2. The constant wind velocity. The wind velocity of 
a hot star normally increases from a value < lkm/s 
in the stellar photosphere to the terminal velocity 
which is reached at a distance of a few stellar radii. 
In polar direction, the B[e] supergiants show a typical 
radiation-driven CAK-type wind with a terminal ve- 
locity of about 2000 km/s. This value was used in our 
calculations. Since the ionization radius in polar direc- 
tions lies at distances of 10 5 — 10 7 AU, the assumption 
of a constant wind velocity, i.e. the terminal velocity in 
these directions, is valid. The only exception is model 
C which is the wind of highest density in our sample. 
Here, the polar ionization radius is reached already at 
a distance of 0.9 AU. For this model a proper inclusion 
of the radial velocity distribution is necessary when the 
correct value of the polar ionization radius is needed. 
In the equatorial direction, in which we are mainly 
interested, our choice of the terminal velocity clearly 
overestimates the real velocity in these regions, since 
we found that the ionization radii here are reached al- 
ready very close to the stellar surface where the wind 
has not yet been accelerated effectively. However, if 
we would include the correct, very slow wind veloc- 
ity, the hydrogen density in equatorial direction which 
is proportional to w _1 would be much larger than in 
our model calculations. The overestimation of the wind 
velocity in our model therefore shifts the ionization ra- 
dius to a larger distance. 

3. The point source approximation for the mean stel- 
lar radiation field. For distances near the star we nor- 
mally would have to take into account effects like limb 
darkening and the size of the stellar disk. Both effects 
lead to a decrease in the mean optical depth at distance 
r from the star. Ignoring these two effects, Eq. re- 
sults in an overestimation of the stellar radiation field 
and in a shift of the ionization radius in equatorial 
direction to a larger distance. 

The assumption of a constant wind velocity that equals 
the terminal velocity and the assumption of the point 
source approximation for the mean stellar raditation field 
result both in a shift of the ionization radius to larger dis- 
tances from the star. Our calculations therefore reflect the 
"worst case" and give an upper limit for the distance of 
the equatorial neutral region. 



7. Comparison with observations from other hot 
stars 

The B[e] supergiants are not the only class of stars that 
shows non-spherically symmetric winds and we mention 
here two galactic stars, the luminous blue variable (LBV) 
AG Carinae and the peculiar B[e] star MWC349A, that 
are also found to have emission from CO and dust coming 
from regions near the stellar surface. 

AG Carinae is a LBV star, even one of the most lumi- 
nous ones with L = 1.7- 1O 6 L (Humphreys et al. HMfy 
and with a variable effective temperature between 9000 K 
and 25 000 K (Lamers HMb Voors et al. 123131) ■ The mass 
of AG Carinae is of order 35 M Q (Vink & de Koter l2UU2}l . 
It is surrounded by a very bright, extended nebula which 
is composed of ionized gas and dust. Most recently, the 
properties of the neutral material around this star was 
subject to a project of Nota et al. I|2()02|l . The authors ob- 
served the lowest rotationally transitions of CO, and with 
their spatially resolved observations they found, most as- 
tonishingly, that the CO emission arises in regions near 
the star. The authors therefore suggested that the CO is 
located in an outflowing disk, just as it is the case for B[e] 
supergiants. Since LBV's have very high mass loss rates 
of order 10~ 5 — 10~ 4 MQyr -1 even at times in between 
individual outbursts (e.g. Leitherer et al. H994|l . the wind 
scenario in AG Carinae might equal those of the B[e] su- 
pergiants, and a neutral disk extending nearly down to the 
stellar surface might surround AG Carinae. Our model cal- 
culations shown above speak in favour of such a scenario. 

The other interesting star is the peculiar B[e] star 
MWC 349 A. Although known for nearly 100 years, its evo- 
lutionary state is still under debate because its spectra 
show characteristics of a pre-main sequence B[e]-type star 
as well as characteristics of a B[e] supergiant (Lamers et al. 
HUH). MWC 349 A has a luminosity of 3 • 10 4 Lq, a mass 
loss rate of 1.2 ■ 10" 5 MQyr" 1 (Cohen et al. ITUgB) . an ef- 
fective temperatue of 35 000 K, and a mass of 25 - 30 Mq 
(Thum et al. 119921 I1994|) . This star is known to possess 
a Keplerian rotating disk. Strong CO first-overtone band 
emission has been observed which must be located near 
the star (Kraus et al. 12000(1 . i.e. much nearer than about 
3 AU which is the dust evaporation radius (Kraus 2000). 

These are only two examples that show how common 
the phenomenon of neutral gas nearby hot stars and lo- 
cated in a shielded, outflowing disk is. 

In our calculations presented in this paper we took an 
empirical density distribution that accounts for the ob- 
served quantities like the density ratio between equator 
an poles, the mass loss rate and velocity in polar direc- 
tions, and the density concentration towards the equator. 
The physical reason for such kind of density distribution 
is not yet fully understood. Since the B[e] supergiants are 
known to be fast rotators, a differential rotation of the 
star with a drastically increase of rotation speed towards 
the equator might lead to the observed angle dependent 
mass flux. There exist several theoretical models to de- 
scribe the concentration or focussing of the material to- 
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wards the equatorial region. The most promising are the 
rotationally induced bi-stability mechanism (see Lamers & 
Pauldrach 1991. Pelupessy et al. 2000) and wind compres- 
sion leading either to a wind compressed disk (Bjorkman 
& Cassinelli EMS! Bjorkman ITM5I Owocki et al. ITSM! 
199(3 H998fl or to a wind compressed zone (e.g. Ignace et 
al. 199G). Unfortunately, none of them is able to explain all 
the observed properties of the polar and equatorial winds 
self-consistently. A better understanding of the physics of 
these outflowing disks and of the ionization structure in 
non-sphcrically symmetric winds is therefore necessary. 



8. Conclusion 

Our main result is that the equatorial winds of B[e] super- 
giants, which are formed due to a high mass flux of low ve- 
locity, can become neutral very near the stellar surface. In 
these neutral disk-shaped regions, which are shielded from 
direct stellar radiation, molecules and dust can form. The 
existence of molecules and dust in the disks around B[e] 
supergiants is already known from observations, but it is 
astonishing that our calculations result in neutral disks 
that extend down to the hot stellar surface. 

These results coincide very nicely with the most re- 
cent spatially resolved observations of strong CO emis- 
sion. This emission seems to originate in regions very near 
the hot LBV star AG Carinae. Its existence can only be 
understood if a neutral disk shields the material from stel- 
lar dissociation. For this star, a scenario as predicted by 
our model calculations, seems to be a natural way to al- 
low CO molecules to exist very close to the hot star. The 
same holds also for the peculiar B[e] star MWC349. This 
star also shows strong CO emission coming from the inner 
parts of a neutral Keplerian rotating disk. 

In addition we find that in polar directions the approx- 
imate ionization radius calculation for winds with an r~ 2 
density distribution turns out to be wrong because one 
uses the physically wrong assumption that there must be 
a turnover from ionized to neutral at a distance smaller 
than infinity. In fact, for large distances, the ionization 
balance equations become independent of r and the ion- 
ization fractions become constant so that no neutral region 
can be found as long as the density in the wind does not 
change, e.g. when reaching the shockfront where the wind 
hits the interstellar medium. 

Although our model calculations are restricted to a 
simple hydrogen plus helium wind, they nevertheless mir- 
ror the general behaviour of the ionization structure in 
non-spherically symmetric winds of hot stars with a high 
mass flux in equatorial direction. Moreover, since the stel- 
lar parameters in our calculations are for a very lumi- 
nous B[e] supergiant, the results are valid for the com- 
plete group of B[e] supergiants, because for less luminous 
stars the number of Lyman continuum photons decreases 
so that the wind becomes neutral even closer to the star. 



Appendix A: Derivation of the classical 
Stromgren sphere 

To show that the classical Stromgren sphere relation fol- 
lows from our more general formulae we start with Eq. 114JI 
for a pure hydrogen nebula of constant density 

oo 

4rrW(r)(l - q{r)) J ^^e~ T "dv = aq(r)Nn (A.l) 

In addition we use the point source approximation, i.e. 
W(r) = l/4:(R^/r) 2 , and multiply both sides of the equa- 
tion by 4nr 

oo 

(1 - q(r)) J ^Rl^j^er^dv = aq{r)N^r 2 (A.2) 

Now, we integrate both sides over the radius and take into 
account the condition for q(r) which is defined by Eq. (|lfi|l . 
The optical depth can then be calculated 



T u =Ja u N H (l-q(r))dr = J a„N R dr = a„N H (r~R s ) (A.3) 

R, Rs 
and Eq. I|A.2|) becomes 

oo oo i?s 

J^Rl^^ Je~ a » N ^ r ~ R ^drdv = aN^ Jdr (A.4) 



hv 

vq Rs 



R, 



The integration over the radius on the left-hand side re- 
sults simply in 

oo 

e -a v N H (r~R,s) dr = _J_ (A. 5) 

Rs 

so that we end up with 



hv a^Nn 3 



or, if we use L v = A-kR 2 -kB v and the assumption Rs S> R* 



^-dv = aN^Rl (A.7) 
hv 3 



which is the classical formula for the Stromgren sphere in 
a constant density nebula. 

For a density distribution of the form iVu(r) = A ■ r~ 2 
the Stromgren sphere is 



which can be derived from our general equations in a sim- 
ilar way. 
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